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CHAOS IN A SPATIAL EPIDEMIC MODEL 1 

By Rick Durrett and Daniel Remenik 

Cornell University 

We investigate an interacting particle system inspired by the 
gypsy moth, whose populations grow until they become sufficiently 
dense so that an epidemic reduces them to a low level. We consider 
this process on a random 3-regular graph and on the d-dimensional 
lattice and torus, with d>2. On the finite graphs with global disper- 
sal or with a dispersal radius that grows with the number of sites, 
we prove convergence to a dynamical system that is chaotic for some 
parameter values. We conjecture that on the infinite lattice with a 
fixed finite dispersal distance, distant parts of the lattice oscillate out 
of phase so there is a unique nontrivial stationary distribution. 

1. Introduction. The inspiration for this paper arose almost 20 years 
ago. The first author had recently moved to Ithaca, New York and the 
Northeastern United States was in the midst of a gypsy moth infestation. 
For all of one summer, he and his wife destroyed egg masses, picked larvae 
off of trees and put bands of sticky tape to keep the larvae from climbing the 
trees. When the next summer came, the outlook for their trees seemed bleak, 
but suddenly all of the larvae were dead or deformed, a victim of the nuclear 
polyhedrosis virus, which spreads through the gypsy moth population once 
it becomes sufficiently dense. 

To model this process we use dynamics that occur in discrete time with 
each site in some graph Gn either occupied or vacant. The number of nodes 
in Gn will be an increasing function of N which tends to infinity. Two 
processes occur alternately: growth and epidemic. 

Growth. Gypsy moths lay dormant in the winter as eggs, so no occupied 
site survives to the next time period but gives birth to a mean (3 > 1 number 
of individuals. Each individual born at x is sent to a site randomly chosen 
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from Nn(%) Q Gn, the growth neighborhood of x, which contains all of the 
nearest neighbors of x in the graph but in general will be larger. 

Epidemic. With a small probability ctN an infection lands at each site. If 
the site x is occupied an infection starts which spreads from x to all of its 
occupied neighbors in the graph and continues until all sites in the connected 
component of occupied sites containing x are wiped out (observe that the 
larger the cluster of occupied sites, the more likely it is to be wiped out 
by the epidemic). It is assumed that the epidemic occurs rapidly so it is 
completed before the next growing season. 

Our goal is to study this process on a random 3-regular graph and on a 
discrete torus of dimension d>2. The second graph is more realistic from 
a biological point of view, but the first one is easier to deal with because 
explicit formulas are available. In both cases, infections will be transmitted 
along edges connecting neighbors. Observe that if we assume that ctN 
then only components with 0(1/ a^) sites will be affected by epidemics. In 
site percolation on an regular tree of degree 3 and on Z d there is phase tran- 
sition from all components small to the existence of an infinite component 
at some density p c . On the random 3-regular graph and the torus this phase 
transition produces one giant component of size 0(n). Thus we expect that 
the density of occupied sites will increase until p> p c , at which point a large 
epidemic occurs and reduces the density to a low level and the cycle begins 
again. We will show that in some cases this leads to chaotic behavior of the 
densities. 

1.1. Mean- field growth on a random 3-regular graph. To work our way 
up to proving results about this system and the corresponding process on the 
torus we begin with the case in which Gn is a random 3-regular graph with 
N nodes, that is, a graph chosen at random from the set of graphs with 
./V vertices all of which have degree 3 (N must be even). We will denote 
this random graph by Rn and we will condition on the event that Rn is 
connected. It is known [see Janson, Luczak and Rucinski (2000)] that the 
probability that Rn is connected tends to 1. We choose this graph, not 
because it reflects reality, but because Rn is locally a tree, so we have 
explicit formulas for the percolation probabilities. To have a simple process 
in which the number of occupied sites at the beginning of the growing season 
is a Markov process, we let A/tv(x) = Rn for all x. As we will see, in the limit 
as iV~ — > oo the result is a very interesting dynamical system. 

To guess what this limiting system must be, observe that if we assume 
that the density of occupied sites before the growth step is p, so the expected 
number of occupied sites is pN, then the expected density after the birth 
step is 
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Now the random 3-regular graph looks locally like a tree in which each vertex 
has degree 3 (we will refer to this tree as the 3-tree) . Proceeding heuristically, 
in the limit N — > oo, each occupied site survives the epidemic if and only 
if it is not in the giant component of the percolation process on the 3-tree 
defined by declaring open the sites that are occupied after the growth step. 
Thus if the density before the epidemic is p, the density gr(p) after the 
epidemic (the T in the subscript is for tree) is exactly the probability that 
the origin is open in this percolation process but it does not percolate. The 
threshold for the existence of a giant component is p c = 1/2, so if p < 1/2 
then grip) =p. 

To compute the density for p > 1/2 we need to compute the percolation 
probability on the 3-tree. Throughout the rest of the paper, whenever we 
say percolation we mean the event that the origin is an infinite cluster of 
occupied sites. We start by noting that for site percolation on the binary 
tree (which is an infinite rooted tree where each vertex has two descendants, 
so all vertices have degree 3 except for the root which has degree 2) the 
percolation probability #bin(p) satisfies 

#bin(p)=p(l-(l-0bin(p)) 2 ) 

since for this event to occur the origin must be occupied and percolation 
must occur from one of the two neighbors. Solving gives 

0bin(p) = ^ = 2--. 

P V 

On the 3-tree the probability of percolation is then 

0T(p)=p(l-(l-0Un(p)f) 

since the site must be occupied and percolation must occur from one of the 
three neighbors. Thus for p £ (1/2, 1] 

(1 \ 3 (1-p) 3 

gr{p) = IP(0 is occupied, Co < oo) = p — 9t{p) =P\ 1 = 5 — • 

VP / P 

Let ao be the solution of 1 — e _/3a ° = 1/2 [i.e., ao = (log 2)//?]. Combining 
the formulas for / and gT we see that the limiting dynamical system should 
be the one defined by the function 

(l-e-PP, 0<p<a , 

Observe that Ht is continuous in [0, 1] . 

We are interested in properties of the iterates of hxip)'- 

• If/9 < 1 then f(p) <p for all p > and thus hj,{p) decreases to as k — > 00. 
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hr(p) 




Fig. 1 . Graph of hr with f3 — 2 log 3. The point c = h T 1 (ao) will play a role in the proof 
of Theorem 1. 

• If > 1 then starting from a small positive p, f k (p) increases to a unique 
fixed point p* . If p* < 1/2 then we never get an epidemic and hj, (p) in- 
creases to the same fixed point. 

• 1/2 is a fixed point when e _/3 / 2 = 1/2, that is, (3 = 2 log 2. When (3 > 2 log 2, 
we let ai = h T (l/2) = e~ 3/3/2 /(l - e _/3 / 2 ) 2 . Eventually the iterates of h T 
lie in the interval [ai,l/2], and once they reach this interval, they stay 
there (see Figure 1). 

Hence if j3 < [3 C = 2 log 2, = f(p) for all p and the epidemic part 

of the dynamics is not seen in the limiting system. If (3 > (3 C then Ht{p) < 
1/2 < f(p) for p> oq. 

Figure 2 shows the orbits of the system as a function of j3. We plot hj^(p) 
for 501 < k < 550 to remove the initial transient. Note that the system pro- 
ceeds directly from a stable fixed point to a "chaotic phase" rather than via 
period doubling bifurcations of the type occurring in the quadratic maps 
rx(l — x). To say in what sense the behavior is chaotic, we will use two 
results of the theory of discrete time dynamical systems. The first result, 
which we include here for convenience, is commonly referred to as "period 
three implies chaos": 

Proposition 1.1 [Theorem 1 in Li and Yorke (1975)]. Let F:J — > J 
be a continuous function on a real interval J and assume that there is point 
a£ J such that 

F 3 (a)<a<F(a) <F 2 (a). 
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Then: 

(a) For every k = 1, 2, . . . there is a point in J of period k, that is, a point 
r G J such that F k (r) = r but F J (r) / r for < j < k. 

(b) There is an uncountable set S C J containing no periodic points such 
that: 

(b.i) For ei>en/ p,q £ S, p^q, 

\i-msvv\F N (p) - F N {q)\ >0 

AT— >oo 

and 

liminf ^(p) -i^fa)! =0. 
TV— >oo 

(b.ii) For every p G S and any periodic point q G J, 

limsup|F%) -F N (q)\ > 0. 

7V-+oo 

We will say that F is chaotic if i* 1 satisfies the conditions (a) and (b) 
above, (b.ii) rules out convergence to periodic orbits, while (b.i) shows that 
all the points in S have different limiting behaviors. 

Theorem 1. 

(a) The dynamical system defined by the function Ht ■ [o>i, 1/2] — > [a±, 1/2] 
is chaotic for every (3 > 2 log 2. 




0,5 I \'y 2 2.5 3 



Fig. 2. Orbits of the system (h^{p))k>o started at p = 0.1. The x-axis has the values of 
(5 used in the simulations, while the y-axis has hij r {p) for k = 501, . . . , 550. 
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(b) If (3 € (2 log 2, 2.48] then the system has an invariant measure, fi = p, o 
Kj} , which is absolutely continuous with respect to the Lebesgue measure. 

Simulations suggest that (b) actually holds for all f3 > 21og(2). 

Now we come back to the process running on Rjy. We will denote our 
process by rj^ , with rj^ (i) = 1 if i is occupied at time k and rj^ (i) = if 
not. The density of occupied sites at time k will be denoted by p^ : 

i=l 

The initial distribution t)q of the process will always be assumed to be a 
product measure with some density p G [0, 1] (so, in particular, p$ converges 
in probability to p) . In the preceding discussion we argued heuristically that 
converges to the deterministic system defined by /it- The next result 
shows that this is indeed the case: 

Theorem 2. Assume that Gn = Rn an d that the infection probability 
of the epidemic satisfies 

on log 2 N > 00. 

AT— >oo 

Then the process (p^)k>o converges in distribution to the (deterministic) 
orbit, starting at p, of the dynamical system associated to hj-- 

The above convergence means that (p^)k>o converges in distribution to 
a deterministic process whose paths are given by the orbits (/i^(p))fc>o- 

1.2. Local growth on the d-dimensional torus. Turning now to a more 
realistic setting, we consider the process running on the d-dimensional torus 
(Z mod N) d , for d > 2, which we will denote by T^r. The case d = 2 is the one 
relevant to gypsy moths, but it is no harder to prove our results in general. 

To add some more realism and make our process more interesting, we will 
take now the growth neighborhoods Mn(x) to be smaller than Tjy- We let 

Nn(x) = {y G T N : < \\y - x\\oo < r N } 

(here the difference y — x is computed modulo N) and take the range r^ to 
be such that r^ — > 00. (We remark that on TV we are considering the L 1 
distance; in particular, two points x,y £ T^r are neighbors if \\x — y\\\ = 1). 

We start as before by guessing what the limiting system should be. To do 
this we will assume for a moment that r^ = 00 for all N, so we are back in 
the case of mean-field growth of the previous subsection. The growth step 
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behaves exactly as before: if p is the density of occupied sites before the 
growth step, then the density after is 



The behavior of the epidemic step in the limit N — > oo is analogous to the 
one in the random 3-regular graph: if p is the density of occupied sites 
before the epidemic, then the density gh{p) after (here the subscript L is for 
lattice) is the probability that the origin is open but does not percolate in 
a site percolation process in 7L d . 

Unlike the case of percolation on the 3-tree, we do not have an explicit 
formula available for the percolation probability in Z d , but we still know 
some qualitative properties. Letting Cq be the percolation cluster containing 
the origin and 



we have that there is a p c G (0, 1) (p c « 0.593 in d = 2) such that 9l(p) = 
for p < Pc, 0l(p) is strictly increasing for p > p c , and 6l(p) is infinitely 
differentiable at every p ^ p c [see Theorem 8.92 of Grimmett (1999)]. We 
also have that 



so is infinitely differentiable at p^p c and gi{p) = P for p < p c . 

As before we let h^p) = 5l(/(p)) and (5 C be the value of f3 solving p c = 
1 - e -/3pc , that is, 



(Pc Ri 1.516 in d = 2). Observe that gi{p) G (0, 1) for p G (0, 1) so, in par- 
ticular, Hl(p) > for p > 0. Our next result holds under an hypothesis on 
the percolation function which might seem strange at a first look, but which 
holds in d = 2 and is expected to also hold in 3 < d < 6. 

Theorem 3. Suppose that 



Then there is an e > such that for every j3 G (/3 C; /3c + £) the dynamical 
system {h k L (p))k>® has an invariant measure which is absolutely continuous 
with respect to the Lebesgue measure. 

We believe (and simulations suggest) that the result holds for all (3 > (3 C - 
As Yuval Peres pointed out to us, it is easy to show that (1.2) holds in d = 2 




L (p)=P(|C o | = oo) 



g L (p) = P(0 < < oo) = P(|Cb| < oo) - P(|C | = 0) =p - 9 L {p) 




(1.2) 



lim 0'l(p) = oo. 

pipe 
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Fig. 3. Density process (pj! (p))k>o running on the two-dimensional torus with mean-field 
growth, depicted for k — 501, . . . , 550 for different values of the parameter (3 (similar to 
Figure 2). Here N = 500. 

using Russo's formula and the fact that the expected number of pivotal sites 
goes to infinity as p [p c in two dimensions. This argument would obviously 
work in other dimensions too if we knew that the expected number of pivotal 
sites blows up at p c . This should be the case in 3 < d < 6 because it is 
expected that 9l{p) ~ C{p — p c ) 7 as p [ p c with 7 < 1 in d < 6, 7 = 1 in 
d > 6, and with logarithmic corrections in d = 6 [see, e.g., Chapter 9 of 
Grimmett (1999)]. 

Our next goal is to show that the process p^? on the torus Tat converges to 
the deterministic orbit of the dynamical system defined by hr, . The processes 
r]£ and p^? are defined in this case exactly as for the random 3-regular graph; 
see (1.1) and the preceding lines. If we consider the case of mean-field growth 
[i.e., Mn(x) = TV] then the result follows from the same arguments as those 
we will use to prove Theorem 2 (the proof is actually simpler because we 
do not have to prove that the torus looks locally like Z d ). Figure 3 shows 
part of the trajectories of p^ in the case of mean-field growth. But, as 
we mentioned, we want to deal with the more general case Mn(x) = {y € 
Tat : < ||y — x||oo < rjy} with r^r — > 00. The result does not seem to be true if 
we do not take — > 00. As Figure 4 shows, the graph of {(p^ , Pk+i)i ^ ^ 0} 
does not correspond to any function. This difficulty dissappears as N — > 00 
if we take — > 00 at an appropriate rate. 

We will assume the following on a at and r^: 
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Fig. 4. Graph of against pjt+i on the two-dimensional torus with N — 750 and 
rjv = 50. The graph clearly does not correspond to a function. 

and 

OiN r N ► OO. 

For instance, we could take = -/V 7 and a^v = N~ s for some < 5 < 7 < 1. 

Theorem 4. Assume that Gn = r ^N, with d>2, and that the number 
of individuals to which each occupied site gives birth to during the growing 
season is a Poisson random variable with mean (3. Then the process (p^)k>o 
converges in distribution to the (deterministic) orbit, starting at p, of the 
dynamical system associated to hi,. 



1.3. Local growth on Z . We now consider the case in which r^ is con- 
stant. Figure 5 shows that when r^ = 5 the fluctuations in the density of 
occupied sites decrease as the system size increases. Figure 6 shows a picture 
of the process running on the torus of size 450 x 450 with rjy = 5. As this 
picture suggests the density stays constant because different parts of the 
lattice oscillate out of phase. 

Theorem 5. Consider the process running in 7L d with d>2. If r^ = L 
and L is sufficiently large then there is a nontrivial stationary distribution. 
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Sketch of the proof. The key to the proof is that the density of 
occupied sites after growth is at most f(l) = 1 — e' 13 so after the epidemic 
there will be a positive density of occupied sites. Let 5 = (1 — e~^)e -4 ^ be 
the probability that a site is occupied and has four vacant neighbors. Divide 
space into squares of side L/2 and declare that the square is occupied if at 
least a fraction 5/2 of the sites are. If L is large enough and T is chosen 
suitably then the set of occupied squares at time nT dominates oriented 
percolation with p close to 1 and the result follows from standard "block 
construction" arguments [for an account of this method see, e.g., Durrett 
(1995)]. By order of the Associate Editor further details are left to the reader. 
□ 

The remainder of the paper is devoted to proofs. The proof of Theorem 
1 is given in Section 2. If you get bored with all of the algebra and calculus 
involved you can skip to Section 3 where the proof of Theorem 2 is given. 
The proof of Theorem 3 given in Section 4 and the more complicated proof 
of Theorem 4 in Section 5 rely on ideas from Sections 2 and 3, but are 
independent of each other. 

The authors would like to thank referee Nicolas Lanchier for his careful 
reading of the paper which resulted in a number of corrections and clarifi- 
cations. 

0,35 

0.3 
0.25 

0.2 
0,15 

0.1 
0,05 


100 200 ,100 400 300 BOO 700 800 

k 

Fig. 5. Sequence of densities pf? of the process running on the two-dimensional torus 
with local interactions in the epidemic step for N — 500 and N = 1500, both with rjv = 5. 
As this graph suggests, the fluctuations of the density process get small as N grows if the 
range rjv is held fixed. 
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Fig. 6. State of the process at time 200 on a torus of size 450 x 450 (black dots are 
occupied). In this simulation, (3 = 2.25, rjv = 5 and the infection probability at each site is 
5 ■ 10 -6 . This picture corresponds to an intermediate state of the process, after an epidemic 
event wiped out a big cluster but the process has had time to grow back. 



2. Proof of Theorem 1. By Proposition 1.1, to obtain (a) it is enough to 
prove that there is a point cG [ai, 1/2] such that 

h T (c) <c<h T (c) <h 2 T {c). 



In our case we can take 



c = f 1 (a ) = llogf 



P °V/?-log2, 

(see Figure 1). Observe that since do < 1/2, c = [3~ l log((l — ao) _1 ) < x 
log 2 = ao • Hence 

M c ) =/( c ) = a o, 

4(c) = /(a ) = i 

and 

h\.{c) = h T {l/2) = a 1 . 

It is clear then that c < hr{c) < h^(c). To see that h^(c) < c we need to 
show that a\ < /~ 1 (ao)> that is, that 

e-W 2 1/ 
< — log 



(1 _ e -/?/2)2 - p b V/3-log2 
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or, equivalently, that 

/ /?e~ 3/3/2 \ 3 

(2.1) UP)= ^ J ^-^)< U i 1) = J JL- 2 

for all (5 > 2 log 2. If you look at the picture of these two functions it seems 
clear that the inequality holds, but the proof is not as simple as the picture 
suggests. We will divide it into two parts. 

First, assume that 3 G (2 log 2, 1.75]. We will show that 

(2-2) M(3 )<A-JL<MP)- 

log 2 

To get the first inequality let 

^-3/3/2 



a{(3) 



(l_ e -/3/2)2- 
A simple calculation gives 

9e /3 _ 4e /3/2 + ! 



a"(0) 



4e 5/3/2 _ 16e 2/3 + 24e 3 P/ 2 - 16e^ + 4e^/ 2 

and we claim that this quotient is positive. Indeed, it is easy to see that 
the numerator is positive, while putting a = e@l 2 the denominator becomes 
4a 5 — 16a 4 + 24a 3 — 16a 2 + 4a, so dividing by 4a we need to show that 

w (a) = a 4 - 4a 3 + 6a 2 - 4a + 1 > 

for all a > 2. Observe that w'(a) = 4a 3 - 12a 2 + 12a - 4, so w'(2) = 4, while 
w"(a) = 12(a - l) 2 > 0, so w'(a) > for all a > 2. Since w(2) = 1 we deduce 
that w(a) > for all a > 2 as required. Hence a is convex, and thus so is 
4>\ = exp(er(-)). Since 

<£i(21og(2)) = 2 = 4- and 0^1.75) « 1.4518 < 4-^—^1.4753, 

Y v &v " log2 Y y } log2 

the convexity of 4>i gives the desired inequality. 
To get the second inequality in (2.2), observe that 

(fe(2 log 2) = 2 = 4- and $,(2 log 2) 



log2 rzv ° ' log2 

Therefore, since this last quantity is exactly the slope of the line appearing 
in the middle term of (2.2) and since <p2 is strictly convex, we deduce that 
(j>2 (P) is larger than this slope for every Q > 2 log 2 and thus the inequality 
holds. 

Now we assume that 3 > 1.75. Using the Taylor expansion of the functions 
1/(1 — x) and e x about x = we get that (2.1) is equivalent to 

v /log2\" > v 1 / 3e-^ 2 
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so it is enough to show that 



log2\ n 1 ( (3e-^l 2 



(3 J ~ n\ \{l-e-P/ 2 f 

for all n > and /3 > 1.75. The inequality holds trivially for n = 0, so by 
induction it is enough to prove that 

log 2 1 /3e" 3/3 / 2 



/3 - n (1 - e-/ 3 / 2 ) 2 
for all n > 1 or, equivalently, for n = 1. That is, we need to show that 

-3/3/2 

for all /3 > 1.75. To see that this holds we observe that the derivative of the 
left side with respect to (3 is 

/3e-^/ 2 (3/?e^ 2 -4e^ 2 -/? + 4) 

2 ( e /3/2 _ I) 3 

We claim that this quotient is negative for (3 > 1.75. Indeed, the denominator 
is clearly positive, so we only need to show that 

W (J3) = 3/Je^ 2 - 4e /3/2 - (3 + 2 > 

for (3> 1.75. This is easy, because w'{(3) = 3e^/ 2 (l + /3/2) - 2e^/ 2 - 1 > e 13 / 2 - 
1 > and w(1.75) ~ 5.28. Thus the left side of the (2.3) is decreasing in (3, 
and then the inequality holds because its value at (3 = 1.75 is approximately 
0.6523 < log 2. This finishes the proof of (a). 

To get (b) it is enough to show by Lasota and Yorke (1973) that 

( 2 - 4 ) pr 1* J(4)'(P)I>1 

p6[oi,l/2]\{ao} 

for [3 G (2 log 2, 2.48]. The idea of the proof is the following. We find an 
explicit formula for (hj*)' and use it to compute numerically its infimum 
on [oi,l/2] \ {ao} for every (3 in a certain grid of (2 log 2, 2.48]. Due to 
monotonicity properties of the derivative of Kt the numerical computation 
of the infimum is exact (up to floating-point numerical errors which are 
small enough for our purposes) for any fixed (3. We then show that (h^) 1 , 
as a function of (3, has a Lipschitz constant that ensures that the infimum 
is larger than 1 for every (3 between subsequent points in the grid. We will 
do this step by step. 

We begin by computing {h\)' . For p £ [oi,ao), h' T (p) = f'(p) = f3e~@ p , 
while for p £ (ao, 1/2], 

TW (l-e-^) 2 (l-e-^) 3 (l-e-^) 31 J 
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This gives an explicit formula for h' T . On the other hand, 

(2.5) (4)'(p) = h' T (h 2 T (p))h' T (h T (p))h' T (p). 

Putting these two formulas together we get an explicit expression for (hj,)' . 

Now observe that h' T is decreasing in [a\,ao) and increasing in (ao,l/2]. 
Indeed, hj,(p) = f"(p) = —(3 2 e~@ p < on the first interval, while on the sec- 
ond one h'^p) = <#(/(p))/'(p) 2 + g' T (f(p))f"(p), so since /' > 0, /" < 0, 



9t(p) 



l-i 

P 



+ 3p[ - 



and 



i 



1 + 



1+1 

p 



1 



1-1 

p 



>o, 



<0 



we get that h'^(p) > for p £ (ao, 1/2]. This means by (2.5) that {hj)' is 
monotone in each interval of constancy of its sign. These intervals are given 
by the partition of [oi,l/2] defined by the preimage of ao under hj,. We 
deduce that 



pe[oi,l/2]\{ao} 



inf min{|(4)'(p-)|,|(^)'(p+)|}, 

p£h T J (a )U{a 1 ,l/2} 

where the superscripts — and + indicate left and right derivatives, respec- 
tively. Using this observation we can compute numerically the infimum in 
(2.4) for any given (3. We did this for every (3 in a grid of width 2 • 10 -6 of 
(2 log 2, 2.48], and we obtained that the infimum is larger than 1.002 at each 
of these values of (3. Figure 7 shows a graph of the values obtained. 

The last step is to make sure that the infimum in (2.4) stays above 1 for 
every (3 E (2 log 2, 2.48]. We will write hr{jp,l3) to indicate the dependence 
of Ht{p) on the value of the parameter (3. Our goal is to find a bound for 
\ g^g p hj^(p, (3)\. Observe that by the product rule and (2.5), if |J^t(p>/3)| < 

Mi and \^h T (p,(3)\ <M 2 for all (3 G (2 log 2, 2.48] and p £ [ai, 1/2] \ {a } 
then 

(2.6) 



d 2 



d(3dp 



h 3 T (p,(3) 



< 3MfM 2 



for all such (3 and p. We already computed \-§^hriPi For p£ [ai,ao), it 
equals (3e~^ p which is smaller than 2.48 for each j3 < 2.48. For p E (ao, 1/2] 
we know that h' T is negative and increasing, so 

e -3/3/2 



0_ 

dp 



h T ( P ,(3) 



< 



d_ 

dp 



-/3/2A3 I 



-3/3 + (3e 



-3-2.48/2 
< —s 4 • 2.48 



1.923. 
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Fig. 7. Infimum of \{h%)'(jp)\ on the relevant interval for (3 £ (21og2,2.6). The compu- 
tation was done for each j3 on a grid of width 2 ■ 10~ 6 on this interval, as explained within 
the proof of Theorem 1. The infimum lies above 1.002 for (3 G (2 log 2, 2.48] . 



Thus if we take Mi = 2.48 the desired inequality holds. Now for p 6 [oi, ao), 
d 2 



dj3dp 
For P e (ao.1/2], 
d 2 

h T ( P ,P) 



kr(p,P) 



df3 [P ' 



|(1 -P 2 )e 



-Pp\ 



< 1. 



dp dp 



-3/3p 



-(- 



-3/3 + Pe 



-p P - 



-Pp 



< 



(i - e-foy 

e ~Pao 

{1-e-P/ 2 ) 



-|(9/3p ~ 3)e" 2/3p + (4 - 4Pp)e- 3/3p + (J3p - l)e" 4/3p | 



-(Upp + 8) < 



1.48/2 y 



-(14-2.48/2 + 8) » 49.73, 



(1-e- 

ifpe (2 log 2, 2.48]. Thus if we take M 2 = 49.73 we get by (2.6) that \ -^h^(p, 
P)\ < 917.6. 

The bound we just obtained implies that for any fixed p £ [a±, 1/2] \ {ao} 
the function p \— > -^hxip, P) is Lipschitz and its Lipschitz constant is at 
most 917.6. Now fix P £ (2 log 2, 2.48] and let P' be the point in the grid of 
(2 log 2, 2.48] on which we computed the infimum in (2.4) which is immediatly 
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before (5. Then for any p £ [a±, 1/2] \ {ciq}, 



d_ 

dp 



h 3 T (p,/3) 



> 



d_ 

dp 



h 3 T (p,(3'] 



> 1.002 - 917.6|/3 - (3'\ > 1.002 - 917.6 • 2 • 10~ 6 w 1.0001. 
This completes the proof of (2.4). 

3. Proof of Theorem 2. To prove this result it will be enough to study 
the one-step transition probabilities for p^? . Recall that in the growth step, 
since here 7Vjv(x) = Gjy, every site becomes occupied with probability 1 — 
(1 — f3/N) pN ~ 1 — e~^ p , where p is the starting density of occupied sites. For 
simplicity we will assume that the occupation probability of each site after 
the growth step is exactly 1 — e~@ p , and then in the proof of the theorem 
we will say how to remove this assumption. 

Abusing notation, we will also let -q^ stand for the set of occupied sites in 
the process. T)k+i/2 wm denote the intermediate state of the process between 
r]fc and rjj^ +1 after the growth part of the dynamics has been run but before 
running the epidemic. We will denote by {0, . .. ,N — 1} the set of nodes of 
Rn. B(i,r) will denote the set of sites in Rn at distance at most r from i 
(here the distance between two points i and j is defined as the number of 
edges in the shortest path going from i to j). 



Let r]i be the set of occupied sites after the epidemic is run on rjy 2 
ignoring infections coming from a distance greater that (log 2 -/V)/5. Define 
rji\/N. Recall that we are assuming that 

aArlog 2 N — > oo. 



Lemma 3.1. 

E(|pf-pf|) >0, 

uniformly in the initial density p. 

PROOF. By translation invariance, and observing that rji(i) > T]i(i) for 
all i £ Rn, 

E(|pf -pf|)<^ E E(|^(i)-r ? f(i)|)=P(0G^)-P(0Gr ? f) 

= P(0 G fjf \ Vi) < (1 - a7v) 1/51og2 N « e~ l ^ aN N —> 0. 

The second inequality above follows from the fact that if is in rj^ but not in 
r)i , then there must be an open path in rjy 2 going from to 8B(0, (log 2 N)/5), 
and all sites in this path must have not been infected. □ 
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Now let 

Hn = {i G Rn '■ B(i, (log 2 N)/5) is a finite 3-tree}. 

By a finite 3-tree we mean a finite tree where all nodes have degree 3 except 
for the leaves which have degree 1. The next lemma says that Rn looks 
locally like a 3-tree: 



Lemma 3.2. 



E( ±-\R N \H N \) =F(0£H N ) >0. 



Proof. A random 3-regular graph is a special case of a graph with 
a fixed degree distribution and can be studied using techniques in Section 
3.2 of Durrett (2007). To explore the subgraph JS(0, (log 2 JV)/5) of R N , let 
Rq = 0, Aq = {0} and Uq = {1, . . . ,N — 1}. These are called the removed, 
active and unexplored sites, respectively. If A n ^ then to go from time n 
to n + 1 we pick a site i n from A n according to some given rule and let 

R n+ 1 = R n U {i n }, 

A n+ l = (A n \ {i n }) U {j G U n :j ~ »}, 



J7 n+ i = [/"„ \ {j G U n :j~i}, 

where j ~ i here denotes that j and i are neighbors. For n < 3iV 1 / 5 /2, |A n | < 
3iV 1 / 5 /2 + 2, so the probability of a collision (i.e., that when we examine the 
neighbors of i n we see a site already in A n ) at some time is at most 

3 U r l 3N 1 / 5 /2 + 2 

2--N 1 ' 5 >0. 

2 N 

Now suppose that when choosing the sites i n we choose those at distance 1 
from first, then those at distance 2, etcetera. Then by time 3-/V 1//5 /2 we 
will have investigated all points within distance (log 2 iV)/5 of 0, and if we 
see no collision, then we will know that the subgraph B(0, (log 2 N)/5) is a 
tree. □ 



Lemma 3.3. Let Co be the cluster containing the origin in a site perco- 
lation process on the 3-tree, and let P p denote the law of this process when 
each site is retained independently with probability p€ [0,1]. Then for any 
k N T oo, 

sup |P p (diam(C ) < oo) - P p (diam(C ) < k N )\ >0. 

pe[0,l] N^co 
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Proof. The result follows from the fact that any increasing sequence 
of continuous functions on [0, 1] which converges pointwise to a continuous 
function on [0,1] actually converges uniformly to that function [see, e.g., 
Theorem 7.13 in Rudin (1976)]. We only need to observe that P p (diam(Cb) < 
oo) and P p (diam(Co) < fcjv) are continuous on [0, 1] as functions of p, and the 
latter is increasing in N and converges pointwise to the former as N — > oo. 
□ 

Lemma 3.4. 

E(hrj?nH N \) >h T (p), 

uniformly in the initial density p. 

PROOF. Observe that since G rj^ implies that G rjy 2 = r]y 2 , 

(3.1) 

= P(0 G 7?f|0 eH N n 7 ? f /2 )P(0 G H N )F(0 G rff /2 ). 

By Lemma 3.2, P(0 G fljv) — > 1 uniformly in p, while by our assumption, 
P(0 G rtf /2 ) = 1 - e-l*. 

For the other term on the right-hand side of (3.1), we only need to look 
at the configuration of t]^ 2 inside -6(0, (log 2 iV)/5), on which, conditional 
on the event {0 G Hn}, the graph looks like a finite 3-tree. Thus, we can 
construct the random variables (rji (®))n>o conditioned on {0 G Hjy H r)^ 2 } 
on a common probability space in the following way. Let T be the set of sites 
in an infinite (rooted) 3-tree and consider a site percolation process on T 
with each site being open, independently, with probability 1 — e~^ p . We will 
call Co the corresponding percolation cluster containing 0. We also consider 
a collection (i3^ v )j g j j 7v>o of independent Bernoulli random variables with 
P(Bf = 1) = oat. With this, the random variable 77^(0), conditional on the 
event {0 G Hjy n f?f/ 2 }, can be constructed as 

~tv (0) f 1, if B? = for all 1 G C n B(0, (log 2 N)/5), 
\ 0, otherwise. 

It is clear that this construction gives the right conditional distribution for 
r?f(0). 

— 1/2 

Now let In = log2( a Ar )• Observe that In < (log 2 N)/5 for large N, so 
we have that 

p(oe^|oeffyn»$ 2 ) 
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= P(0 G rj? ,diam(C ) < l N \0 £H N n r$ 2 ) 

(3.2) 

+ P(0 G r/f , Zat < diam(Co) < 5 log 2 AT|0 G #jv n r?^ 2 ) 
+ P(0 G r/f , diam(Co) > i log 2 A^|0 eH N n rff /2 ). 
For the first probability on the right-hand side we have that 
P(0 G 7/f ,diam(C ) < Zjv|0 G #tv n r/^ 2 ) 
< P(0 < diam(Co) < l N \0 £H N D rj^ /2 ) 

g T (l - e-to) 



< diam(Co) < 00 10 is open) 



This convergence is uniform in p thanks to Lemma 3.3. On the other hand, 
since any subset of T with diameter n has at most 1 + 3 • 2 n_1 < 3 • 2 n nodes, 
we get that 

P(0 G rff , diam(Co) < l N \0 G Hn n r]? /2 ) 

= F(Bf = Vi G C , diam(Co) < l N \0 £H N n rff f2 ) 

= E((l - aiv) |Co1 , < diam(Co) < /jv|0 £H N n rff /2 ) 

> (1 - a N ) 3a " 1/2 F(0 < diam(Co) < l N \0 G Hn H r^ 2 ) 

; gr(l-e-^) 
~~ ' 1 - e"^ 

, -1/2 

uniformly in p by the same reason as above and because (1 — ayv) N ~ 
e -3^ _, 1. We deduce that 

P(0 G rf , diam(Co) < Zjv|0 G i?jv H < /2 ) — ^f," ^ , 

uniformly in p. For the second probability on the right-hand side of (3.2) we 
have that, since P(0 eH N D r]^ /2 ) > C = (1 - e~ /3p )/2 for large enough N, 

G rj?,l N < diam(Co) < ± log 2 iV|0 G Fjv n V1/2) 
< C -1 P(Ziv < diam(<7 ) < ±log 2 iV) 
= C" 1 [P(diam(C ) > l N ) - P(diam(C ) = 00)] 

- C" 1 [P(diam(C ) > § log 2 iV) - P(diam(C ) = 00)] 
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uniformly in p, again by Lemma 3.3. For the last probability in (3.2) we 
simply observe that 



P(0 G rj?, diam(Co) > ± log 2 N\0 € H N n r^ /2 ) 



< (1 - a N 



.l/51og 2 JV 



-1/5ojv log 2 N 



0. 



The previous calculations and (3.2) imply that 



N 



1 - e-fo ' 

uniformly in p. Putting this together with (3.1) we get the result. □ 



Proof of Theorem 2. By Karr (1975), it is enough to prove that 
Po P an d that given any sequence pn in [0, 1] converging to some p' 6 [0, 1] , 
the sequence p± , with t)q started at a product measure of density pjy, 
converges weakly (or, equivalently, in probability) to 1it(p')- 

The first part is straightforward. For the second part we will assume, for 
simplicity, that pn = p' for all N and, moreover, that each site is occupied 
with probability 1 — e~@ p after the growing season. The general case follows 
from the facts that 1 — (1 — (3/N) p converges uniformly as N — > oo to 1 — 
e —pp for p' £ [0, 1] and that, by the preceding lemmas, all the convergences 
we will prove below are uniform on the initial density p. 

Observe that by Markov's inequality, given any e > 

F(\pf-h T (p')\>e)<-E(\p^-h T (p')\), 



so 



(3.3) 



Ml0l>e)<-E(|pf 
+ -E 

£ 



N 



1 

+ -E 

£ 



~N 
Pi 



?j?nH N \ 



N 



l 

+ - 

£ 



E[-\^nH N 



h T (p') 



Lemmas 3.1 and 3.4 imply that the first and last terms on the right-hand 
side of the inequality go to as N — ► oo. The second one also goes to since, 
using Lemma 3.2, 

1 



E 



Pi 7 



±\v?nH N \ 



<E(-\R N \H 



N 



N 



0. 



To deal with the third term, observe that 
Var(l^fn^l) 
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N-1N-1 

i=0 j=Q 

< eH N x H N :B(i,(log 2 N)/b)nB(j, (log 2 iV)/5) + 0}| 

= |{(t, j) €H N xH N : < 2(log 2 iV)/5}| < 2iV • N 2 '\ 

Hence, by Jensen's inequality, 



E 



±\?j?nH N \-v(±\rj?nH N 



2 <Var(^l|^nif w 
2N ■ N 2 / 5 

^— >0 - 

We deduce from (3.3) that p± converges in probability to /jt(p'). ^ 

4. Proof of Theorem 3. As in the case of the 3-tree, we let ao be the 
solution of /(ao) = p c [i.e., a = log(l/(l - p c ))/P] and ai = /il(p c ) (see 
Figure 1 for a sketch of these values in the case of the 3-tree) . It is enough 
to prove, by Lasota and Yorke (1973), that there is a K G N such that 

(4-1) r mf |(/#)'(p)|>l. 

Fix any f3\ > f3 c . Since a\ is bounded away from for (3 G (/3 C5 /3i), there 
is a K 6 N such that minjfc £ N: f k (a\) > p c } < K — 1 for any such /3. In 
particular, since ao is always less than p c we deduce that given any j3 G 
(PaPi) and any p G [ai,p c ], the -ff -tuple (p, /iz,(p), ■ ■ • _1 (p)) contains at 
least one point in (ao,p c ]- 

Now recall that /" < 0, so /' attains its minimum on the interval [ai,ao] 
at ao, and at this point its value is (3(1 — p c ). Thus for every [3 G (/3 c ,/3i), 
this minimum is larger than (3 C (1 — p c ). Since gi{p) = p for p G [ai,ao] we 
deduce that 

\h>L(p)\ > PcO- ~Pc) for all p G [ai,a ]. 
Now using the fact that ao T Pc as J. /3 C , we can choose given any e > a 
/% G (Pc,Pi) so that /(p c ) -p c = /(p c ) - f(a ) < e for any (3 G ((3 C ,(3 2 ). Since 
(1.2) implies that 

0l(p) = 1-0l(p) ►-<*>, 

PlPc 

we can choose a small enough e, so that 

\h' L (p)\ = \g' L (f(pMf(p)\>^{mi-Pc)}- {K - 1 \l} 

for all p G (ao,p c ], and thus this inequality holds for all f3 G (/3 C , /in- 
putting the previous arguments together with the fact that 

(fcf )'(p) = ^(fcf- 1 (p))^(fcf- 2 Cp)) • • • h' L (p) 
we deduce that (4.1) holds for all (3 G (/3 C) /?2)- 
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5. Proof of Theorem 4. Given i G TV and m G N we will write 

B(i, m) = {j G Tjv : ||^ — j||oo < m } and V(m) = (2m + l) d = \B(i,m)\ 

(here and in what follows all differences i — j for i,j G HV are computed 
modulo N). Define, for k G N, 

lb-*l|oo<r-]v 

and 

G^(e) = {*€T JV :|d^(*)-/ l i(p)|< e }. 

(i) is the density of occupied sites in the growth neighborhood of i, while 
G 1 ^ (e) can be thought of as the set of "good sites at time k" where a site 
is said to be good at time k if the density of occupied sites in its growth 
neighborhood at that time is close to the desired value h\(p). The proof of 
Theorem 4 will depend on the following proposition: 

PROPOSITION 5.1. Fix ei,£2 > and k G N and assume that 
(5-1) ±E(\T N \G»(5 1 )\)<6 2 . 



Then if 5± and 5 2 are small enough and N is large enough, 

l 1 E(\T N \G^ 1 (e 1 )\)<e 2 . 



This result will allow us to give an inductive proof of Theorem 4. We will 
need thus the following lemma: 



Lemma 5.2. Given any 5 > 0, 



Proof. By translation invariance, 

E(\T N \ G$(5)\) = £ P(t * G»(5)) = N d F(\d$(0) -p\>5). 

Since ^(d^ (0)) = p, Chebyshev's inequality and the fact that (by definition) 
V(r N )d£(0) is the sum of V^r/v) independent Bernoulli random variables 
with success probability p imply that 

P(|<(0) -p\ > S) < 62v l {rN)2 V(r N )p{l - p), 
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so 



1 , ^jv«- - 



^ 3 E(|T^\Gg'(*)|)<^ 5 ^ P (l-p)— 0. D 

Now we turn to the proof of Proposition 5.1. Many parts in the argument 
will be similar to those in the proof of Theorem 2 and the lemmas that 
preceeded it, so we will skip some details. We begin with some preliminary 
results. Throughout this part, and until the proof of Theorem 4, we fix 
k, Si, 62, £1, £2 and assume that (5.1) holds. 

Observe that since each occupied site i sends a Poisson[/3] number of 
births during the growing season, each to a site chosen randomly from 
B(i,rN), we can equivalently think of each occupied site i as sending a 
Poisson[/3/F(r7v)] number of births to each of its V(rj^) neighbors at a dis- 
tance smaller than r^. Hence during the growing season, each site % receives 
Y.\\3-i\\ 00 <r N v lk U) Y j,i births, where (Y i:j ) i:j& j N are i.i.d. Poisson[/J/^(r;v)] 
random variables. Conditional on d^{i), this last sum is distributed as a 
Poisson[(fj^ random variable. We deduce that we can regard the grow- 
ing season as taking place as follows: 

Given r)™ , each i will be in ri^ +1 , 2 with probability equal to the probability 
that a Poisson[d^ (i)0\ random variable is positive, that is, with probability 

The Poisson random variables above are taken to be independent of each 
other. 

Let In = \Aiv /&n and observe that 

and oinIn = ^/ajjrjq — ► 00. 



r N ^Ja N r N 

-1/2 



We let rj^ +l be the configuration obtained from rj^ +1 ^ 2 by ignoring infections 
coming from a distance greater than 

Lemma 5.3. 

-L E^+iW-dWD— ^ 0. 



In particular, 



N d - TV-,00 



E(|pf +1 -pf +1 |)^0. 



Proof. By translation invariance, and repeating the arguments of the 
proof of Lemma 3.1, we get that 

^ E e(|^ +1 (o - Ci (01) = jip E p (*Ai« * 
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= !P(OeCi\'Jfc+i)<(l-«Jv) is 

w e" Qjv/jv — ► 0. □ 

Before continuing, it is useful to give an explicit construction of the ran- 
dom variable rjj^ +l (0). Consider a collection X = (X(i)) ieIi d of i.i.d. random 
variables with uniform distribution in [0, 1] and, given rj^ , construct f)^ +1 ^ 2 
as follows: 

1 &l/2(*) = l ;r( . )>e -<(i)- 

Observe that with this choice, ^(Vk+in^f) = 1) = 1 — e -|9d fc W as required. 
We will call Cq the open cluster in T]^ +1 , 2 containing 0. Define (B^) ieZ d^ >0 
as in Section 3 and set 




!' if ^4-1/2 (°) = 1 and B i = fOT a11 * G C B ( ^ l N), 

0, otherwise. 



This construction gives the right distribution for ?/^ 1 (0). 

We introduce another modification of rjj^ +1 : let fjj^ +1 be the configuration 
obtained from rj^ in the same way as Vk+i' exce Pt that in the growing season 
we replace ilk+1/2 by the configuration fj^ +1 , 2 defined by 

d/2« = i x( . )>e -^ ( p) 

(using the same family of variables X). That is, Tjk+l/2 corresponds to run- 
ning the growth step as if the density of occupied sites in the ball of radius 
tn around each site was exactly h k L (p). p^ will denote the density of oc- 
cupied sites in this modified process, that is, p^ = \r}J?\/N d . We will call 
Co the open cluster containing in the site percolation process in all of Z rf 
constructed from the collection of random variables X with each site being 
open with probability 1 — e~^ h ^ p \ 

Lemma 5.4. Given any e > 0, if 5\ and 62 are small enough, then 

E(|pf +1 -pf +1 |)< £ . 

Proof. The idea behind the proof of this result is the following. By 
(5.1), the density of occupied sites is close to around most sites. If 

this holds for some site i, then in a box around i the density must still be 
close to this. We then prove the result by comparing rjj? +1 and ffj^ +1 with 
processes in which the outcome of the growth step is replaced by product 
measures of sligthly smaller and slightly larger densities. 
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To get started we observe that 

<4 E E(iCi«-C+i»i)=p(d(o)^d(o)) 



(5.2) 



N d 

iGTjv 



< P(Ci(0) / Ci(0),0 G + P(0 i G?{6i)) 

< P(C+i(0) / C+i(0), G G^(fc)) + S 2 , 

where in last bound we used (5.1). To deal with the last probability we first 
observe that given any i G 5(0, Zjv), 



^) = ^3 E ^(i) 



^(0) + ^ E ^(j) 

^ N> j£B(i,r N )\B(0,r N ) 

E 



V ( r N) jeB (0,r N )\B(i,r N ) 

N \B(i,r N )\B(0,r N )\ 

- dk (0)+ VW) 

and thus, since the cardinality in the last term is largest when i is at any of 
the 2 d corners of the hypercube B(0,In), we have that for some C > 

df(o)|<c^«^o. 

V(r N ) r N 

We deduce that 

n^ +1 (0)^Ci(0)>0eGf(^)) 

< P(Ci(0) ± Ci(0), |4« - >*£(p)l < 25i V* G 5(0, Zjv)> G G^(<5i)) 
+ P(|d^(i)-/i|0)| >2«Ji for some iG 5(0, l N ),0 e G% (S^) 

< P(Ci(0) / Ci(0), |4« - &£(p)l < 2<5i Vi G 5(0, Zjv)) 
+ P(|df (i)-d£ (0)| ><5i for some iG 5(0, Zjv)) 

+ P(|^(0)-/ l i(p)|>5 1 ,0GG^(5 1 )) 
= P(Ci(0) / Ci(0), |^(») " h k L (p)\ < 25 l Vi G 5(0, l N )) 

for large enough N. 

Next, we introduce the following notation: £', 2 will be the set of open 

sites in a site percolation process in Z d with each site being open with 
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probability 1 — e~^ q for q G [0, 1] constructed from the family of random 
variables X. In other words, we put £jy 2 (*) = ^-X(i)>e-P* f° r each i G 7L d . We 

also let C T^r be the configuration obtained after running the epidemic 
step on £jy 2 n HV (this is done on the torus HV, so we take into account 
the periodic boundary conditions of the torus while running the epidemic), 
using the variables (B^)i & j N , and ignoring infections coming from a distance 

greater than Zjv- Observe that with these definitions, Vk+1/2 = £1/2^ 

and 77^ ) _ 1 = £^ L ^ ,N . The key fact is the following: 

P(Ci(0) / Vk+M, \d%(i) ~ h k L (p)\ < 25 1 Vt G 5(0, to)) 

<p(e^ (p)+2<5liiV (o)=o,^ (p) ' JV (o) = i) 

(5.3) + p(^)-^(o) = o,^; 2 w (o) = 1) 

+ P(^«-^(0) = 1,^^(0) = 0) 

+p(4 w+25i (o)=i,4 w - 25i (o)=o). 

To see that this is true observe that (i) — h\{jp)\ < 25\ for all i G 5(0, to) 
implies that 

1 _ e -/3(fc*»-25i) < x _ e -/3<W < X _ e -0(fc*(p)+«i) 

for all i G 5(0, to), and thus 

4 (P) ^ 251 n 5(0, to) c n 5(0, J*) c 4 (p)+2<51 n 5(0, to). 

Assuming this, we have that rj^ +1 (0) = and rj^ +1 (0) = 1 implies that 
£^'^(0) = ^ 2 {p) (0) = 1, and either ^ +1/2 (0) = 0, which implies that 
Zu^~ 2Sl (0) = 0, or rjg +1/2 (0) = 1 but there is an infection in n 5(0, to), 
which implies that £^^ +25l,JV (Q) = 0. Similarly, rj^ +1 (0) = 1 and rj^ + i = 
implies that £*y 2 P +2<51 (0) = 1, ^ L ^' N — o, and there is no infection in 

n 5(0, to), and thus ^^"^'^(O) = 1 whenever ^j^" 2 * 1 (0) = 1. 
To finish the proof we need to bound the probabilities on the right-hand 

side of (5.3). boi the hrst one, since £y 2 Q qy 2 — £1/2 ' we na ve 

that if denotes the size of the cluster containing in the configuration 
given by £, then 

p(^ (p)+25l,7V (o) = o,^ (p),7V (o) = i) 
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< P(^ (P)+251 '>) = 0, # (P) '>) = 1, #£f +2Sl < oo) 

+ P(^ (p) '>) = 1, = oo) + P(#4 (P) < #4 (P)+251 = oo). 



The first probability on the right-hand side is bounded by 
(5.4) 



F^ liP)+25l ' N (0) = 0, 4 (P)+251 (0) = 1, M^ P)+251 < oo) 



<E(1-(1-M*4 (P ^ 

which goes to by the dominated convergence theorem. The second one 
goes to as well because it is bounded by (1 — aA?) /jv ~ e _QjvZjv . The third 
one equals 

e L (h k L (p) + 26 1 )-Q L (h k L (p)), 

which is less than e/2 for small enough 5± by the (uniform) continuity of 
the percolation probability 0l(p) for p G [0, 1]. The other two probabilities 
on the right-hand side of (5.3) can be bounded similarly, yielding 

P(Ci(0)^w(0).0eGf(^))< £ 

for large enough N and small enough <5i. Putting this together with (5.2) 
gives the result. □ 

Lemma 5.5. 

l%ti)-^ +1 (p)Ho. 

Proof. This proof is similar to that of Lemma 3.4. First we observe 
that 

E(^ +1 ) = P(o g Cilo g Ci/ 2 )iP(o e C1/2) 

(5.5) 

= P(0Gr?f +1 |0Gf?f +1/2 )[l-e-^W] 

and 

(5.6) P(0G7?f +1 ,diam(C )=oo|0Gd /2 )< (I - a N ) lN & e~ aNlN — >0. 
Now 

P(0 G »?f +1 ,diam(Co) < oo|0 G rj% +1/2 ) 

(5.7) =P(0 G d,diam(C ) < l N \0 G f?f +1/2 ) 

+ P(0 G rj£ +1 ,l N < diam(Co) < oo|0 G C+1/2) 
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and, trivially, 
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(5.8) 



P(0 G In < diam(Co) < oo|0 G 7?f +1/2 ) 

< P(/jv < diam(Co) < oo|0 G C1-1/2) — ► 0. 
On the other hand, 

P(0 G 7?f + l> diam (^o) < MO € C1/2) 

= P(Bf = Vi G C n B(0,/ Ar ),diam(C ) < Zjv|0 is open) 
= E((l - a 7V ) |ConB(0 '^ )l ,diam(C ) < l N \0 is open) 
= P(diam(C ) < l N \0 is open) 

- E(l - (1 - a iV )l ConB ( '^)l,diam(C7o) < l N \0 is open). 
The second expectation is positive and bounded from above by 
E(l - (1 - a N )\ C °\, |C | < oo|0 is open), 

so it goes to as N — > 00 by the dominated convergence theorem as in (5.4). 
Thus 



lim P(0 G ^ +1 ,diam(C ) < l N \0 G fj% +1/2 



/v 



/V 



P(diam(Co) < oo|0 is open) 



_ P(0 < diam(Co) < 00) _ g L (l - e~^W) 
~ 1 _ e -Ph k L ( P ) ~~ 1 _ e -Ph k L { P ) 

Putting this together with (5.7) and (5.8) we get that 

h k L +1 (p) 



P(0 G 7^ +1 ,diam(C ) < oo|0 G f% +1/2 ) 

and thus by (5.5) and (5.6) we obtain 

\E(p% +1 )-h k L +1 (p)\- 

as required. □ 



I_ e -Ph k h (p) 



Proof of Proposition 5.1. 



1 



^E(|T^ \ G» +1 ( £l )\) = P(0 i G£ +1 ( £l )) = P(|4Vi(0) - hf(p)\ > ei ) 



iN 



N 



<-E(i^ +1 (o)-^ + >)D- 
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Hence 

^E(|T^\G^ +1 ( £l )|) 
(5.9) < -[E(|^ +1 (0) -^ +1 (0)|) +E(|^ +1 (0) -E(pf +1 )|) 

+ E(|E(pf +1 ) - + E(|^+i - 

where df +1 (0) = F ^ Ey^^ CfiO')- 

For fixed e > we want to show that each of the expectations on the 
right-hand side of the last inequality can be bounded by e if N is large 
enough and di and 62 are small enough. The bound for the last one follows 
directly from the triangle inequality and Lemmas 5.4 and 5.5. 

For the first one we have by translation invariance that 

E(K +1 (0)-dT +1 (0)|)<-^ £ E(|^ +1 (j)-dO')l) 

V[rN) Hillock 



£ E(|r ? f +1 (j)-^ +1 ( J )|)<s 



for large enough N by Lemma 5.3. 

For the second one we first observe that, again by translation invariance, 
E(^ +1 (0))=E(pf +1 ). Hence 



(5.10) 



E(|^ + i(0)-E(^ +1 )|) 2 
<Var(^ +1 (0)) 



Cov(C +1 «,^ + i(i)) 



V M i,jeB(0,r N ) 

^ 777^2 Km G B ( ^) ■ II* " illoo < In}\ ~ ^4 — * °- 

The bound for the third expectation on the right-hand side of (5.9) follows 
from the exact same argument as previous one. We deduce that 

1 — - . 4e 



— EdT^VG^^)!)^ 
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for large enough N, and thus choosing s < £i£2/4 gives the result. □ 

Proof of Theorem 4. Since [0, 1] is compact, it is enough to prove 
the convergence of the finite-dimensional distributions of p^f , and since our 
limit is deterministic, we only need to prove that 

(5.11) n \ p N_ h k L(p)l>£) Q 

for every k > and e > 0. Proceeding as in the proof of Theorem 2 we have 
that 

n\f$ - h h L (p)\ > e) < -MP% ~Pk\) + -M\P% ~ Pk\) 

(5.12) 

+ ±E(\p»-E(p»)\) + -JE(p%)-h L (p)\. 

By Lemmas 5.3, 5.4 and 5.5, given any v > there are constants S^ 1 , S%~ 1 > 
such that 

(5.13) ^E(|T W \ GtiiSt 1 )]) < 4^ 

implies that the first, second and last terms on the right-hand side of (5.12) 
are each bounded by ve for large enough N. The third term is also less than 
ve for large N, which follows from repeating again the argument in (5.10). 
We deduce that 

(5.14) n\Pk-h k L (jp)\>e)<Av 

for large enough N provided that (5.13) holds. 

Similarly, Proposition 5.1 implies that (5.13) will hold provided that 

M^V) W firm \ n N (xk-2\\\ ^ rfe-2 



N d 



-E(|Tjv \ G ! jt_ 2 (^i )|) < #j 



for some 5\ 2 ,6% 2 > 0. Repeating this procedure inductively we deduce 
that (5.14) holds provided that 

^E(|T,v\GM)|)<<5 2 

for some small 5®, 5® > 0, which holds for large enough ./V by Lemma 5.2, 
and thus (5.11) follows. □ 
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